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Abstract 

We investigate chiral properties of the domain- wall fermion (DWF) system by 
using the four-dimensional hermitian Wilson-Dirac operator. We first derive 
a formula which connects a chiral symmetry breaking term in the five dimen- 
sional DWF Ward-Takahashi identity with the four dimensional Wilson-Dirac 
operator, and simplify the formula in terms of only the eigenvalues of the oper- 
ator, using an ansatz for the form of the eigenvectors. For a given distribution 
of the eigenvalues, we then discuss the behavior of the chiral symmetry break- 
ing term as a function of the fifth dimensional length. We finally argue the 
chiral property of the DWF formulation in the limit of the infinite fifth di- 
mensional length, in connection with spectra of the hermitian Wilson-Dirac 
operator in the infinite volume limit as well as in the finite volume. 
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Typeset using REVTpjX 



I. INTRODUCTION 



A suitable definition of the chiral symmetry has been a long standing problem in lattice 
field theories. Recently an ultimate solution to this problem seems to appear in the form of 
the Ginsberg- Wilson relation |],|2[]. Two explicit examples of the lattice fermion operators 
which satisfies the Ginsberg- Wilson relation have been found so far: One is the perfect 
lattice Dirac operator constructed via the renormalization group transformation |3|,[|] and 
the other is the overlap Dirac(OD) operator []5]|| derived from the overlap formalism J7| 
or from the domain-wall fermion(DWF) |^-[T0[| in the limit of the infinite length of the 5th 
dimension. Since the explicit form is simpler for the latter, a lot of numerical investigations 
18 1 as well as analytic considerations []T9| , pO| have been carried out for the domain-wall 



fermion or the overlap Dirac fermion. 

Recent numerical investigations for the domain- wall fermion [|14],|]1|], however, bring puz- 
zling results, which are summarized as follows. 

Analytic considerations suggest that the overlap or domain-wall fermion works well at 
sufficiently weak gauge coupling or equivalently for sufficiently smooth gauge configurations 
]TP| , |2l][ ] . Initial numerical investigations supported this result [|TT|-|i"5[. On the other hand, 
further investigations indicate that the domain-wall fermion at stronger coupling ceases to 
describe the massless fermion even in the N 5 — ► oo limit [H,I5|, where N 5 is the number 
of sites in the 5th dimension. Analytic results f2T|-p3f in the strong coupling limit are 



controversial. However the latest one |24| also suggests that the DWF does not work in the 
limit. 



It has been argued |14| that this result may be understood by the relation between the 
phase structure of the lattice QCD with the 4 dimensional Wilson fermion [^J and zero 
eigenvalues of the 4 dimensional Wilson-Dirac operator. It is well-known that the zero 
eigenvalues cause a trouble for the domain- wall fermion (or the overlap) [T^|I^J2^] , therefore 
it is natural to consider that the success/failure of DWF depends on the absence/presence 
of the zero eigenvalues. On the other hand, the zero eigenvalues of the 4 dimensional Wilson 
Dirac operator, denoted as D w , is related to the parity-flavor breaking order parameter, 
(gz7 5 r 3 g), where r 3 is the 2x2 flavor matrix. Introducing the external source H coupled 
to qi^r^q, the following relation is easily derived. 



qi^r q 



lim lim — Tr- 

H ^0+ V-^oo V D w + IJ 5 T 3 H 



— lim lim — tr 

H-+0+ V->oo V 



«75 



275 



—i lim lim — tr 
—i lim lim — 

if->0+ V^oo V 



-i lim / d Xp Hw (X) 

H — h) ' J 



D w + il5H D w - i-f 5 H 
1 1 



H w 



iH H w - iH 
1 1 



A n + iH X v 
1 1 



iH 



.X + iH X-iH 
-i I d Xp Hw {X){-2iri)6{X) = -2np Hw (0) 



(1-1) 



where X n and |A n ) are the eigenvalue and the eigenstate of the hermitian Wilson-Dirac 
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operator Hw = 75-tV and Ph w W is the density of eigenvalues of Hw, defined by 



Ph w {>) = Jim — J2 6 ( X ~ X r 



;i.2) 



The expected phase structure of lattice QCD with the Wilson fermion [2S| is given in 
Fig. |l], where g is the gauge coupling and M is the mass parameter. In the region B, 
(gz75r 3 g) 7^ 0, thus the density of the zero eigenvalues is nonzero. According to this phase 
structure, the DWF is successful between < M < 2 in the weak coupling limit of QCD. 
This allowed region of M agrees with the analytic consideration of DWF ||. Once the gauge 
coupling becomes nonzero, the allowed region of M shrinks and moves to larger values. 
This property has been predicted by the mean-field analysis |26| and has numerically been 
observed |14j]. If the coupling becomes large so that (3 = 6/g 2 < /3 C , the allowed range of M 
disappears and the massless fermion ceases to exist in the domain- wall QCD (DWQCD). 

Numerical results mentioned before |l^|,|15j that DWF does not work in the strong cou- 
pling region seem to agree with the above expectation. This interpretation, however, has 
difficulties. Let us summarize the recent numerical result [13], where the quenched DWQCD 
has been investigated with the RG improved gauge action as well as the ordinary plaquette 
gauge action. With the former gauge action the quenched DWQCD works at f3 — 2.6 and 
fails at (3 — 2.2, suggesting that 2.2 < (3 C < 2.6, while with the latter gauge action the 
quenched DWQCD fails even at (3 = 6.0, indicating (3 C > 6.0. One problem is that the latter 
condition [3 C > 6.0 with the plaquette gauge contradicts the previous numerical investigation 
27j, which concludes that the region without parity-flavor breaking (: the allowed region in 
the present case) exists (3 = 6.0. A more serious problem is that the numerical analysis for 



the density of eigenvalues indicates non-zero value of Ph w (0) at any (3 |28| with the plaquette 
action, and the similar conclusion is obtained by the same analysis with the RG improved 
action [p9j] . If this is true, the phase structure in Fig.|I|is incorrect in the quenched QCD: the 
gap of the parity-flavor breaking phase never shows up in the weak coupling region, leading 
to the conclusion that DWF does not work at all in the quenched QCD. 

This complicated situation is summarized as follows. While being consistent with the 
numerical result of the DWQCD for the plaquette action |L4,|15], the numerical analysis for 
the density of the eigenvalues p8[ , |2"9l contradicts the phase structure of QCD with Wilson 
fermion [p5| , pTf and the numerical result of DWQCD with the RG action ||14|| . Furthermore 
the numerical result of the DWQCD with the plaquette action seems apparently inconsistent 
with the numerical result for the phase structure of the Wilson fermion with the same gauge 
action. 

In this paper, we try to resolve the above mutual inconsistency, by analysing eigenvalues 
of the 4 dimensional hermitian Wilson-Dirac operator. In Sec. || the action and the axial 
Ward-Takahashi identity are given for the latter use. In Sec. [ill], we derive the formula for 
the explicit chiral symmetry breaking term of DWQCD, m^, q , in terms of eigenvalues and 
eigenstates of the modified hermitian Wilson-Dirac operator. In Sec. [IV|, assuming the form 
of the eigenstates, the formula derived in Sec. [TTT] is reduced to a much simpler expression 
in terms of the eigenvalues only. The dependence of m^ q on N 5 is discussed in Sec. |V|, for 
a given model of the density of the eigenvalues. In Sec. [VTJ, using the simplified formula, 
we discuss the chiral properties of DWQCD first in the finite volume, and consider it in 
the infinite volume limit. Some comments are also given on the phase structure and the 
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exceptional configurations in the lattice QCD with the Wilson fermion. Our conclusion is 
presented in Sec. [VII . 



II. ACTION AND AXIAL WARD-TAKAHASHI IDENTITY 



We employ Shamir's domain- wall fermion action @,|T(|. Flipping the sign of the Wilson 
term and the domain wall height M, we write 



Sdw = « 4 H VK^j s ) 

x,y,s,s' 

D w 
D 5 



D 



w 



M 
a 



a 



V f V b 
2 v n v n 



(2.1) 
(2.2) 
(2.3) 



where x, y are four- dimensional space-time coordinates, and s, s' are fifth-dimensional or 
"flavor" indices, bounded as 1 < s, s' < N*, with the free boundary condition at both ends 
(we assume iV 5 to be even). The domain- wall height M should be taken as < M < 2 at 
tree level to give massless fermion modes. V{ and V;, are forward and backward derivative 

PA* 

in four dimensions 



1 



x,y 



yj ' 



and V5 and Vg are those in the fifth dimension with the boundary condition 



(2.4) 
(2.5) 
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s+l, 



_ l<s<N 5 ) 
a 5 \ a 5 m f S s/<1 - 8 S>S > (s = N 5 ) 

(s = l) 
(K s < N s ) 



1 

a 5 



d. 



8 S _ S > — 5. 



a 5 m f 8, 



s'N 5 



8-1,8' 



(2.6) 
(2.7) 



where m/ is the mass for the quark field. The light fermion mode is extracted by the 
4-dimensional quark field defined on two boundaries of the fifth dimension, 



q(x) = P L il>(x,l)+P R il>(x,N 5 ), 
q(x) = ^{x, N 5 )P L + ^{x, 1)P R , 



(2.8) 



where Pr/l is the projection matrix Pr/l = (1 ± 7s)/2. The quark mass term induces a 
coupling between these boundary fields through the bare quark mass m/. 

In this paper we investigate chiral property of the domain-wall fermion through the axial 
Ward-Takahashi (WT) identity defined in Ref. [110] for the non-singlet axial transformation 



5 a A ij(x, s) = ie(N 5 + 1 - 2s)TV(x, s), 
Sfflix, s) = -ie(N 5 + 1 - 2s)T>(x, s), 



(2.9) 
(2.10) 
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where e(x) is a sign function of x. The WT identity for some operator O is written as 
follows. 

V b »Al{x)0) = 2m f (P a (x)0) + 2 (J% q {x)0) - (5 a A 0) , (2.11) 

where an axial vector current A a , a pseudo scalar density P a , and an explicit chiral symmetry 
breaking term J 5q at finite N 5 are given by 

A°(x) = £ e(iV 5 + 1 - 2s)-(>(x, S )T a (l - 7/1 )[/ /i (x)^(x + /x, s) 

s=l 2 

-^(x + /i, s)(l + 7m )^(x)T>(x, s)), (2.12) 

P a (x) = g(a;)75T a g(x), (2.13) 

1 /— - N*. 

J 5 » = - i^(x, -f + l)T a P R ^(x } -f) - j,{x, -f)T a P L ^ -f + l) 

= ^'(3)rV(4). (2.14) 

O5 z z 

Here we define the following fields according to [[]] 

iV* iV* iV* 

-y) = PflVfo y) + ^(x, -y + 1), (2.15) 

^(s, ^) = ^(x, ^ + 1)P« + ?(x, ^)P L . (2.16) 
In this paper we consider the identity with O = P b {y) 

K A » pb (y)) = 2m f (p a (*)P b (y)) + 2 (WP b (y) 

-^*,v(q(y){T a ,T b }q(y)) (2.17) 
and measure the chiral symmetry breaking effect by 

""'-fe^iPftfiPior (2 - 18) 



which we call an 'anomalous quark mass' [[14]]. Please notice that we have omitted flavor 
indices since the flavor factors are canceled in the definition of m§ q . 



III. EFFECTIVE THEORY OF DOMAIN- WALL FERMIONS 

We now derive the effective theory of the DWF system by integrating out the heavy bulk 
modes according to Ref . || . Our aim is to rewrite the numerator and denominator of the 
anomalous quark mass ( |2.18| ) in terms of four dimensional quantities and hence to relate m 5q 
to an hermitian Wilson-Dirac operator in the four dimensional theory. Introducing source 
fields to q(x), q(x), if)'(x, N5/2), ip (x,N 5 /2), the propagators necessary for our purpose are 
derived as |J 
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(q(x)q(y)) 



a 5 1 



a 5 1 — am f 



D N l - a 



a 5 



a 5 2 cosh ^a 5 # 



a 5 
a 



2 cosh ^a 5 # 



75, 



(3.1) 
(3.2) 

(3.3) 



where if is a DWF Hamiltonian in the 5th direction, with which the transfer matrix in fifth 
direction is given by 



T = e 



-a 5 H 



(3.4) 



Dn 5 is a truncated Ginsparg- Wilson (GW) Dirac operator which satisfies GW relation in 
N 5 —>■ oo limit. H and D^ 5 are related to the four dimensional hermitian Wilson-Dirac 
operator H w as, 



j\r 5 ~ 

1 + arrif) + (1 — am/)7s tanh — a$H 

2 



~ 1 . 1 + aH' 

H = — log . 

as 1 — aH' 



H' = H\ 



w : 



2 + aj 5 H w ' 
H w = 75 (d w - — "j 



(3.5) 
(3.6) 
(3.7) 
(3.8) 



Here we adopted Borigi's notation [3D| for H. 

The numerator of m^ q can be written in terms of the above propagators, 

X(t) = £(J 5g (x,t)P(f,0)> 

X 

1 

- \ t r -r(r'i.r. f . • ^ 
1 



«5 



£ tj 



v. ( ' ''(. t, ^)q(y, 0)^ 75 ^g(j7, o)^'(f, t, -y) 



x,a,a j3,b 



2 cosh ^a 5 # 



D 



At 



A'.-, 



2 cosh -^a 5 H 



(3.9) 



where J = (af, £, a, a), J = (y,0, fi,b). a, (3 are spinor indices and a, 6 are color indices. \I) is 
an eigen-ket in the coordinate, spinor and color spaces. In the last equality we use a relation 
that 75(-D7Vg)75 = {Djffi . The denominator of m^ q is given by 

Y(t) = ^2(P(x,t)P(y,0)) 

= - £ tr [75 (q(?, t)q(y, 0)) 75 (q(y, 0)q(x, t))] 

X 

al 1 

Tl 1. I ^ 



a 10 (1 — arrif) 2 



££('IK 



— a 



J J 



a 10 (1 — am/) 2 



e e (/ ^; 

x,a,a f3,b 



J) J 



75 

A. 



OJ 75 



(3.10) 
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Note that (I\ J) = for t ^ 0. 

Hereafter we set the bare quark mass rrif = without loss of generality. 



IV. EXPANSION IN TERMS OF EIGENSTATES 

We expand X(t) and Y(t) in terms of eigenstates of the hermitian operator H'. Since 
the numerator ( |3.9| ) has a suppression factor 1/(2 cosh ^-a§H) only small eigenvalues of H, 
or equivalently those of H', contribute to m 5? at large N$. In this case we can expand H' 
perturbatively in terms of the four dimensional hermitian Wilson-Dirac operator Hw and 
can derive the formula which connect m^ q and the eigenvalue of Hw at large N 5 . 

The expansion of X(t) and Y{t) with these eigenstates is done by inserting a complete 
set of eigenstates 

i = 5»h, (4- 1 ) 

n 

where |n) is n-th eigenstate of H' 

H' \n) = X' n \n) . (4.2) 

With this substitution X(t) and Y(t) are written in terms of eigenvalue and eigenfunction. 

In order to further simplify these Green functions as a function of eigenvalue only, we 
employ the following assumptions and approximations for typical form of the eigenfunctions. 
The recent numerical analyses indicate that the eigenvalues of Hw seem to be classified into 
two groups JTP|]^T|; one is a group of small isolated eigenvalues and the other is a group 
of almost continuous eigenvalues above the isolated ones. The eigenvectors associated with 
the isolated eigenvalues are exponentially localized at some center in the coordinate space 
19| , |3TH , and those for the continuous eigenvalues are rapidly oscillating plane- wave functions 



in the coordinate space f31j . This property of eigenvalues is expected to hold also for H', 
since H' can be expanded in terms of Hw perturbatively for small eigenvalues. ^From this 
consideration we assume that the eigenvector space of H' is divided into two subspaces S 
and S, where S is spanned by localized eigenvectors and S is expanded with plane wave 
functions. 



A. Localized eigenvectors 

We adopt two types of the approximation for typical form of the localized eigenvectors; 
completely localized ones and partially localized ones. The completely localized eigenvec- 
tor means that the eigenfunction ip n of the n-th eigenvalue has non-zero value at a single 
parameter (x n , a n , a n ) 

if> n (x, a, a) = 5 XjX J atan 5 atan , (4.3) 

which gives a unit vector in (x, a, a) space. ip n (x, a, a) is given by using the n-th eigenstate 
\n) of an eigenvalue 
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ip n (x, a, a) = (x, a, a\n) , 
£» = A», 
H\n)=X(X' n )\n), 

We assume that eigenvectors with different eigenvalues reside at different points 

(n\x, a, a) (x, a, a\m) =0 for n ^ m. 
In this case S is spanned by a set of basis vectors 









1 




W 



(4.4) 
(4.5) 
(4.6) 

(4.7) 



(4.8) 



(4.9) 



On the other hand, the partially localized eigenvector is non-zero at small range of 
(x,a,a) space. We call this small range of volume h n as S n and S = U n S n . S n is spanned 
by a set of eigenvectors ip^\ where % runs as % — 1, • • • , h n and is given by 



V'n l '(x, a, a) = (x, a, a\n, i) , 
H'\n,i) = \J ] \n,i) . 



We assume that S n with different n has no overlap each other; 

(m,i\x, a, a) (x, a, a\n,j) = 



(4.10) 
(4.11) 



(4.12) 



for m 7^ n for any In the limit of h n — * 1 the partially localized eigenstate is reduced to 
the completely localized one. 

Correspondence between the eigenvalues and eigenvectors of H' and Hw is given pertur- 
batively for small eigenvalues. We expand H' in terms of H w 



2H> = H w + H w ±(-iy( a ^ 
k=i v 1 

and the eigenvalue Ajf and eigenstate \n w ) of H w is defined by 

Hw \nw) = ^ \ n w) ■ 



(4.13) 



(4.14) 



The standard perturbation theory gives the relation between A n , \n) and A„ , \nw) as follows. 



1 



= - - (Af) (n w \ l5 \n w } + 0(X 3 ), 



In) = \n w ) - 



H W 75 \n w ) + 0(X 3 ), 



(4.15) 
(4.16) 



where <p n is a projection operator onto the space perpendicular to the n-th eigenstate: 

(j) n = 1 - \n w ) (n w \ . (4.17) 
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B. Plane- wave eigenvector 



We approximate the plane-wave function with that of the free theory but relax the 
condition for possible momenta p in the finite box. The operator H' of the free theory in 
momentum space is given by 

H\p) = ^75 (*E7^» + > ( 4 - 18 ) 

MP) = s , + {2 %_ Mr (4-19) 

B(p) = ^ + (2^-M)(2 + 2^-M) 
W s 2 + (2 + 2s 2 - M) 2 ' 1 Uj 

where s M = sin (ap^), s 2 = X^ sin 2 (ap M /2) and s 2 = X^ The eigenstate of free if' and if 
becomes 

tf\p,3,a) = \'(p)\p,s,a), (4.21) 

#b,s,a) = A(p)|p,s,a>, (4.22) 

A'(p) = ±-JA\p) + B^p), (4.23) 



A(p) = llogl±^, (4.24) 

where the spinor index s runs as s = 1, 2, 3, 4, and the eigenstate is degenerate in the color 
index a. The eigenvector in (x, a, a) space can be written as 

(x, a, a\p, s, b) = -^=U a {p, s)e ipx 5 a , b (4.25) 

within the subspace S, where U a (p, s) is a normalized eigenfunction of free H'(p), given in 
Appendix [A], and N p is the total number of the plane-wave eigenvectors. The free truncated 
overlap Dirac operator and its inverse in momentum space are given by 



(k, s, a \D Ns \p,t,b) = — (z(7 M ) s ,*C M (p) + 8 s ,tE(p)) S k ,pS a ,b, (4.26) 



DnI 



p,t,b) =2a g2 + 0k,p5 a ,b, (4.27) 



where 



= ^M tanh T as|X(p)l ' (428) 
= ( x + ^M tanh T a5|X(p)l ) • (429) 

In addition we also assume that the overlap of two eigenfunctions from different groups 
vanishes: 
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(n\x, a, a) (x, a, a\p, s, b) = 0. (4.30) 

The eigenvalues of the hermitian Wilson-Dirac operator Hw is also given as a function 
of for the free theory. The free Hw in the momentum space is given by 

H w = 75 (iJM + 2s 2 - M) . (4.31) 

The eigenvalue Xwip) °f f ree Hw is given by 

H w \X w ip)) = V(p) \X w ip)) , (4-32) 



\ w (p) = ±\j s 2 + (2s 2 - M)\ (4.33) 

where \\wip)) is the corresponding eigenstate. Although the eigenstates of X'(p) and Xwip) 
are different, X'(p) of H' is given in terms of Xwip) an d the momentum p as 

X'ip) = XW{P) (4.34) 

IX w {p) 2 + A{2s 2 - M) +4 



C. Simplification of the formula of m^ q 

We first consider a system with iVj completely localized eigenvectors and plane wave 
functions with iVp degrees of freedom for the momentum p. Ni and N p satisfy a relation 
Ni + 4N c N p = 4N C N X , where N x is a number of total sites of our lattice N x = n x n y n z nt- 
Since the total eigenvector space is separated into two subspaces S and S, the complete set 
of this system is given by 



Ni N p 4 N c 

1 = J2\ n ) ( n \ + J2J2J2 b' s ' a ) (p> s ' a l 

n=l p=l s=l a=l 



where |n) G <S and |p, s, a) G <S. 

We expand the numerator of with the complete set in the above. 



= X,(t)+X c (t), 



2 cosh ^a 5 # 



J J 



2 cosh ^a 5 # 



(4.35) 



(4.36) 



where Xi(t) and X c (t) are the contributions from the localized eigenvectors and the plane 
wave function, respectively. After a little algebra, the detail of which is given in Appendix 
M, we obtain 



1 1 Nl 

Xl{f) = ^ r(t) 4Ay^ n ? 1 \2coshf a 5 X(X' n 



(4.37) 



for the contribution from the localized eigenvectors, where Y(t) is the denominator of m^ q . 
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Similarly, taking t — > oo limit, we have for the contribution of the continuous eigenvalues 

2 

F(p), (4.38) 



1 



a» x ^iV n ^ \2coshf a 5 A(p) 



"p p 



where n' is a volume of <S in xyz space, and F(p) is defined as 



F(p) 



4N r 



N p \C 2 (p) + E 2 (p)J ' 
We next evaluate the denominator of m^ q in the similar manner: 



(4.39) 



x,a,a f3,b 



J) J 



where 



N, 



AN C N X 



K 1 )' 



no 



(4.40) 



and 



t— >oo 



(4.41) 



(4.42) 



These results imply 

Y(t -> oo) 



WcNs-Nia 10 xyz N p ^ [P> - 



(4.43) 



We can directly see that the function -F(p) is almost constant for all the range of p as 
in Fig. except some discrete points of p. Therefore it is a reasonable approximation to 
assume that F(p) is independent of p. Combining these formula and approximations the 
anomalous quark mass becomes 



m 5q 



lim^opX^) 
lim^oo Y(t) 

1 1 



E 



a 5 4N c N x IV \ v 2cosh^a 5 A(A; 
where p n is a number of degeneracy in A n , 

Pn = E h n ,x( P )- 

p,ct,a 



+ Ep^ 



2 cosh ^-G^Ar 



(4.44) 



(4.45) 



If the degeneracy of p n in the free case is resolved by the presence of gauge fields, we are 
able to reconstruct m 5q by simply summing up all eigenvalues of H', 
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1 1 \ - 

m5g a 5 4N C N X V V 2 cosh f a 5 A(A' n 



(4.46) 



The above result can be generalized to the case with the partially localized eigenstates. 
The formula becomes 

1 1 



a^N c N x ^coshf a^)) + ? P " ( 2 coshf a 5 A n ) J ' ^ 

where the n-th set of localized eigenvectors have the (local) support with the dimensionless 
volume of h n , and the h$ (1 < h$ < 4 4 /i n ) is the enhancement factor for the set, which 
depends on the shape of eigenvectors. See the detail of the derivation in Appendix |C[ When 
the degeneracy in the continuous eigenvalues is resolved with gauge fields we have 




m ^ = t^t I 2^M — — „ r-— I + x 



4N C N X l local n I 2cosh^a 5 A(A;)/ continuous I 2 cosh ^a 5 A(A^ 




(4.48) 



V. MODEL FOR EIGENVALUE DENSITY 

In this section we consider the relation between the asymptotic behavior of m§ q in N5 
and the distribution of the continuous eigenvalues. For simplicity we neglect the effect of 
the localized eigenvalues here and will discuss their effect in the next section. In this simple 
situation we can write m 5g as an integral in continuous eigenvalues 

m ^=L. dXp{x) {Mfx) ' (5 - 1} 

where A is a dimensionless eigenvalue A = a^X, A m j n and A max are minimum and maximum 
of eigenvalues < A m i n < A < A max . Without loss of generality, we consider non-negative A 
only, by taking p(A) = p„(A) + p n (— A), since cosh ^A is the even function. We adopt the 
following three types of p(A) 

p(A) = sjR 2 - (X-R-S) 2 (5.2) 

with 

1. 5 < 0, A min = 0, A max = 2R- \S\, p(0) ^ 0, 

2. 6 = 0, A min = 0, A max = 2R, p(0) = 0, 
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3. 6 > 0, A min = 6, A max = 2R + 5, p(0) = 0. 

A typical form of p(A) for each case is given in Fig. ^| 

A main support of the integral at large N 5 resides near A m i n and asymptotic behavior 
of m 5q is evaluated by expanding p(A) around A m i n and adopting the leading term. Since 
contribution from the larger eigenvalues is negligible we set the integral range to [A min , oo]. 
The integral in A is easily calculated with the following two formulas, 

^ruT""**, (5.3) 

osn— A) — x 

d\e~ x \ a = { , , (5-4) 




which is valid for x > 0, a > — 1. 

For S < the asymptotic behavior of the anomalous quark mass becomes 

m ^iu ay (5,, 



'5 



For 5 = 

m 5g 



where ( is the Riemann's zeta function. For 5 > we have 



We can see that 5 = gives power law even if the zero mode density vanishes p(0) = 0. 
We need a gap at A = in the continuous eigenvalues to realize exponential decay. Note 
that the similar analysis is also made in Ref . f32| . Typical form of m^ q is given in Fig. f| for 
R = 4 and 5 = -0.5, 5 = 0, 5 = 0.5. 



VI. EFFECT OF LOCALIZED EIGENSTATES AND THE INFINITE VOLUME 

LIMIT OF THE SYSTEM 

In this section, using the formula for m 5(J in terms of the eigenvalues, we propose our 
interpretation, which resolves the inconsistency among the numerical simulations for the 
domain-wall QCD, mentioned in the introduction. 

We consider the situation that the density of the continuous eigenvalues has a gap at 
zero: p(A) = for |A| < 5 with 5 > 0. On the other hand, we do not put such a restriction 
on the localized eigenvalues, so that they can become almost zero. In this situation, the 
contribution of the continuous eigenvalues to m 5(? vanishes exponentially in N 5 . In the 
following subsections we will consider the contribution from the localized eigenstates to m^ q . 
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A. N5 — > oo at finite volume 

In this subsection we discuss the behavior of m^ q in the large N5 at finite volume. This 
situation often corresponds to the one encountered in the numerical simulations, and in the 
infinite N 5 limit the domain-wall fermion becomes the overlap Dirac operator, which satisfies 
the Ginsparg- Wilson relation. 

In the finite volume, it is almost impossible to have an exact zero eigenvalue. More 
precisely a probability to have the exact zero eigenvalue is zero, since no symmetry assures 
the existence of it. Indeed no exact zero is numerically found in the evaluations of eigenvalues 
of Hy/. On the other hand, almost zero eigenvalues, |A| ~ 10~ 2 or less, appear, and they 
become smaller as the volume increases. Therefore it may be reasonable to assume that the 
average of the smallest eigenvalue vanishes in some power of 1/V: (|A m i n |) = c$V~ Cl with 
c ,ci > 0. 

In this situation m^ q vanishes exponentially as N 5 increase: 

m590C 4iv^ exp[ ~l^r ] ((U) 

as N 5 — > oo. This means that DWQCD always works in the finite volume: m 5(J vanishes 
exponentially in N&. In other words, as long as the continuous eigenvalues have a gap around 
zero, DWF in the N$ — > oo limit, or equivalently, the overlap Dirac fermion, works well to 
describe the chiral modes in the finite volume, where the smallest eigenvalue of the localized 
modes is small but non-zero. 

According to this consideration, we speculate how DWQCD behaves as the coupling 
constant varies. In the strong coupling, even the continuous eigenvalues have no gap such 
that |A min | = 0, and therefore m 5g vanishes only in some power of I/N5. DWQCD does 
not work in the strong coupling region. Once the gap in the continuous eigenvalues opens 
(|A m in| > 5 > 0) at P > P c in the weak coupling region, m 5q vanishes exponentially in N 5 
and DWQCD in the finite volume works well to describe the chiral symmetry. 

We think that p c < 6.0 for the plaquette action and p c < 2.6 for the RG improved action. 
At first sight this seems to contradict with the numerical data of m^q, which does not show 



the exponential decay in N$, for the plaquette action |14]jl5|| . This contradicted behavior is 
explained as follows. In the intermediate values of N 5 , the continuous eigenvalues give the 
main contribution of m 5(? , so that 

m 5q ~C coat exp[-N 5 \\™£\}, (6.2) 

while as N$ further increases the localized eigenvalues dominates, so that 

m5 9 ~aocaiexp[-A 5 |A^ n al |], (6.3) 

where C cont or C\ oca i is proportional to the number of the modes near A^?* or A^gf, respec- 
tively. At P — 6.0 for the plaquette action, the transition between the former exponential 
behavior with [A ^fff l and the latter with (A^ 1 ! can be seen for A5 = 10 ~ 40 [14,15]. The 
latter behavior looks almost constant in A^ 5 since Aj^ 1 is very small. At P — 2.6 for the RG 
improved action, on the other hand, only the former exponential behavior can be detected 



for N 5 = 10 ~ 24 fl4j| . This suggests that the ratio Ci oca i/C cont is smaller for the RG action 



14 



than for the plaquette action. Indeed it has been numerically found |29] that the number 



of the localized modes near Aj^ 1 is much less for the RG action, while the number of the 
continuous modes near A™^ is similar in the two actions. (However lA^g^l seems larger for 
the RG action.) We then give two predictions, which should be checked in order to test this 
interpretation: the exponential decay, exp[— -/VslA^ |], can be seen at larger N 5 > 60 for 
the plaquette action, and the transition between the former and the latter can be seen at 
7V 5 = 40 ~ 60 for the RG action. 

The distribution of eigenvalues of H w suggests that even (3 = 5.65 for the plaquette 
action or f3 = 2.2 for the RG action is already in the weak coupling region (J3 > /3 C ). Indeed 
it seems that the transition between two exponentials has been observed in the behavior of 
m 5q at these (3 [|4 



B. Infinite volume limit at finite N§ 

In the previous subsection, we argue that m^ q decay exponentially in N$ in the finite 
volume as long as the distribution of the continuous eigenvalues has a gap around zero. The 
exponential decay rate |A^'|, however, may vanish in the infinite volume limit. Therefore, 
it may be the case that m 5q does not vanish exponentially if the infinite volume limit is 
taken before — > oo. Moreover, the value of N$ necessary for the suppression of the chiral 
symmetry breaking, m^g, may increase as the volume becomes larger. This may be a disaster 
from the practical point of view. In this subsection, we will discuss whether the effect of the 
localized eigenstates to m 5q vanishes or not in the infinite volume limit. 

In the case of the hermitian Wilson-Dirac operator, the plane-wave eigenvalues in the 
finite volume become the continuous spectra in the infinite volume. On the other hand, we 
can not predict the nature of the distribution for the localized eigenvalues in the infinite 
volume limit. Although the localized modes are always discrete in the finite volume, it can 
be continuous in the limit as will be discussed later. 

We first consider the case where the number of the localized modes around zero does not 
grow linear in the volume. More precisely, the number of modes which satisfy — e < A < e 
with e > 0, denoted by N(—e, e), is bounded by coeV Cl with Co > and c\ < 1. In this case 
the contribution to m^ q in the infinite volume limit becomes 

focalized _ nm 1 V- e -|A„|7V 5 < N(-6,e) + Q( ^-eN^ 

n:localizcd 

< lim coeV' 1 - 1 + 0(e- eN5 ) = 0(e" eiV5 ), (6.4) 

V^oo 

where the sum is taken only for the localized modes. This result shows that the contribution 
from the localized near-zero modes vanishes in the infinite volume limit and the remaining 
contribution also vanishes exponentially in iV~5, as long as the number of the localized near- 
zero modes is bounded as iV(— e, e) = coeV Cl with c\ < 1. DWQCD works well in this case 
even in the infinite volume limit. 

If the number of the localized near zero-modes is proportional to the volume; N(—e, e) = 
CqsV, the contribution from them remains non-zero in the infinite volume limit: 

1 

V 



focalized Hm V- e -|A„[JV 5 = ^ + Q ^-eN s y ^ 



n:localized 
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DWQCD does not work in this case in the infinite volume limit. 

If N(—e, e) = c eV, the density of state p(A) is (almost) continuous and nonzero at A = 
in the infinite volume limit: p(0) 7^ 0. In other words, a sum of the infinitely many 5 functions 
from the localized modes, normalized by the volume, becomes the continuous function in 
the infinite volume limit. To be more precise, one should first define the integrated density 
of states by 

k(X) = lim -N(-oo,X). (6.6) 

V-+00 y 

Since k(X) is a monotonically increasing function in A, its derivative, , becomes a 

dX 

well-defined measure. The previous statement is equivalent to that pi OC aiized(A) = — ^— is 

dX 

continuous and non-zero at A = 0. 

Let us explain the case that piocaiized(O) 7^ more concretely by using a model of the 
eigenvalue distribution. Suppose that the discrete eigenvalues are uniformly distributed in 
the interval that — 1 < A < 1 and the number of modes is equal to 2cV, so that the average 
interval between two successive eigenvalues becomes l/(cV). In the infinite volume limit, 
we see that the eigenvalues are discrete but dense in the interval. If we calculate m 5q using 
this distribution, we have 

m ^ = Jim tt £ e '~ 5 = J im 77 " = tt ! - e~ N * . 6.7 

Therefore m 5q does not vanish exponentially in N$ if the infinite volume limit is first taken. 

The appearance of the continuous density of states from the localized modes in the 
infinite volume is often observed if a kind of randomness is introduced in the interaction 



Physically the Anderson localization ]3B[ is one of such examples. In QCD one can 
have infinitely many localized modes if infinitely many pairs of instanton and anti-instanton 
with a fixed topological charge exist. In lattice QCD, we have, in addition, very localized 
modes associated with the dislocations. Since the dislocations is local, the number of the 



discrete modes associated with them can be proportional to the volume ]37 , 34 



Although the contribution of the localized modes to m^ q vanishes in the finite volume, 
whether it vanishes in the infinite volume limit depends on the density of states of the lo- 
calized modes in the infinite volume limit. The previous numerical investigations on p(0) 
28| , |29[| suggests p(0) is small but non-zero at V/3 7^ 00 in quenched QCD. Further investi- 



gations, however, are necessary in particular for the volume dependence of p(0), to have a 
definite conclusion on this point whether p(0) is non-zero or not. If p(0) 7^ at V/3 7^ 00, 
it is interesting and important to derive the (3 dependence of p(0) theoretically. Note 
however that p(0) should vanish in the continuum limit for < M < 2, since no zero mode 
appears in the free theory. 



C. Phase structure with the Wilson quarks 

In this subsection we consider the relation between the above interpretation for the 
behavior of m^ q in DWQCD and the parity-flavor breaking in the Wilson quark action in 
the quenched approximation. 
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We first show that the contribution of the localized modes to the parity-flavor breaking 
order parameter vanishes in the infinite volume limit, in the case that N(—e, e) < c^eV C1 
with Ci < 1. In this case we have 

(^i7 5 r 3 ^) localized = lim lim \- £ ~ 2H 



mlocauzed n 



< lim lim — 

H^OV^oo V 



-2H _ -2H 



H 2 ,ft. e 2 + H 2 



2 i 2H 

< lim lim coeV Cl_1 — + lim lim — C{e)V- 



h^o v^oo H h^o v^oo V e 2 + H 2 

where I_ n; |A n |> e 1 = C(e)V. The first term vanishes in the V —>■ oo limit, while the second 
one vanishes in the H — > limit. 

We next show that the contribution of the localized modes to the parity-flavor breaking 
order parameter remains non-zero in the infinite volume limit, if the density of such states 
is non-zero at A = in the infinite volume limit: N(—e, e) = c^tV . To see this, we again use 
the uniform distribution of eigenvalues between —1 and 1, used in the previous subsection. 

9 cV o it cV ffrV 

(^7 5 r 3 ^) localized = lim lim - V — — = -4c lim lim V 



H_.oy-.ooy ^ (^) 2 + # 2 M^» n fi« 2 + (M) 2 

_ 7rcoth(7ri7cV) 1 



—4c lim lim 

fl->0 V-kx> 



2#cV 



+ A 



-2c7T, (6.9) 



where lim#_+ limy^oo A < lim^^o ttH/ (2a/ 1 + H 2 ) = 0. 

The above two considerations show that the correspondence between the failure of 
DWQCD in the infinite volume limit and non-zero order parameter of the parity-flavor 
breaking in the quenched Wilson fermion still holds even in the case that the localized 
modes dominate in near-zero eigenvalues. The previous investigations p8| , ^9| suggests that 
the gap in the phase structure in Fig. [I] closes even at (3 > (3 C in the quenched QCD with 
the Wilson quark. In this case we expect that non-zero value of pi OC aiized(0) in the region A 
is much smaller than p C ontinuous(0) in the region B, the true parity-flavor breaking phase |2S 



so that p(0)iocaiized may be too small to be detected by measuring the pion mass. This may 
be a reason why the previous investigation indicates that the region A exists at (5 = 6.0 for 
the plaquette action ^7j. Finally it is noted that the region A with piocaiizcd(O) = should 
always exist in full QCD, since appearances of the localized near-zero modes are suppressed 
by the fermion determinant, det Dw = det Hw- 



D. Exceptional configurations 

Finally let us comment on the exceptional configurations appeared in the quenched QCD 
with the Wilson-type quark action. It is sometimes observed at small quark masses in the 
quenched simulations that the pion propagator receives the anomalously large contribution 
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from a particular configuration, which is called the exceptional configuration and is removed 
from the statistical average. In ref. [^] the several configurations, which give the W shape 
in the pion correlation function in t space, have been observed in the region A, where the 
gap (absence) of the party-flavor breaking phase is expected. It has been very difficult to 
understand the W shape propagator, since it means that the correlation increase as the 
separation t increases. Now we argue that this may be understood by the localized modes. 
The pion correlation function is defined by 



{if)ij 5 r a ip(x) • i)i^T a i)(y)) = (x 



1 






l 












x^ 



£Wn)f <A n |y)5>|A,)i (A,|x) 



i A n Az 



(6.10) 



where 0„(x) = (x\X n ) etc. If some localized eigenstate has a very small eigenvalue A n on 
some configuration, the contribution of the mode becomes very large: 



(i/ji-f 5 T a i)(x) ■ i)i^T a i)(y) = — <pn{x)<p n {x) ] (t)n{y)(t)n{y) ] + other contributions. (6.11) 

K 

This is the exceptional configuration. Furthermore suppose the eigenstate 4> n (x, t) is localized 
at x = (x n ,t n ). Summing over x and y, and taking x = t and y = 0, the time correlation 
is dominated by 



c (t) = E 72 ' *)^(^ tfMv, o)0n(£ o) f , 

- - A„ 

x,y n 



(6.12) 



and C(t) has a peak around t = t n . Together with an enhance factor 1/A^ for small 
eigenvalues this explains the W shape behavior of the pion propagator. 



VII. CONCLUSIONS 

In this paper we first derive the formula for the chiral symmetry breaking term m^ q in 
DWQCD in terms of the (modified) hermitian Wilson-Dirac operator in 4 dimensions. Using 
several simplifications and approximations we explicitly write down the formula of m^ q in 
terms the eigenvalues only. 

The important observation is that there are two different types for the eigenvalues of the 
the hermitian Wilson-Dirac operator. One is the continuous one, which corresponds to the 
plane-wave in the free theory, the other is the localized one, associated with the instanton 
or dislocations. We argue that the effect of the latter one to the physical observable such 
as m 5q or the parity-flavor breaking order parameter vanishes in the infinite volume limit, 
unless the number of near zero localized modes increases linearly in the volume. On the 
other hand, the effect remains non-zero but small in the infinite volume limit, if it linearly 
increases. 

The message of this paper is that domain-wall fermion or overlap Dirac fermion should 
work well to describe the vector-like chiral symmetry at weak coupling, (3 > j3 c , where the 
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gap in the continuous spectra opens, as long as the volume is finite. The small eigenvalues 
appeared in the numerical simulation seem to belong to the localized one, hence their con- 
tribution vanishes in the N$ — > oo limit in the finite volume. Whether domain- wall/over lap 
fermions are successful or not in the infinite volume, however, depends crucially on the 
distribution of the localized eigenvalues in the limit. 

In practice one can make such small eigenvalues appeared in the finite volume larger 
by hand, without changing physical observables [^,^]. This is a little costly. If the effect 
of chiral symmetry breaking is small enough or no dependence of observables on N 5 is 
detected, one may instead perform simulations at large but numerically affordable value of 
N 5 . Except the quantities very sensitive to the small eigenvalues such as m 5g , N 5 dependences 
are expected to be rather week in general. This expectation is indeed true in the case of the 



quantum hall effect HO 
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APPENDIX A: PLANE WAVE 



The eigenvector ( |4.25| ) of the free Hamiltonian ( f4.18|) is given by combining positive and 
negative energy eigenvector u a (p, r) and v a {p, r) 



(H') a pVp(p,T) = -\X'(p)\v a (p,r), 
{H') aP up{p,T) = \X'(p)\u a (p,r), 



(Al) 
(A2) 



where a, f3 are spinor indexes and r runs r = 1,2. We adopt the following combination in 
this paper; 

U a (p, 1) = v a (p, 1), U a (p, 2) = v a (p, 2), U a (p, 3) = u a (p, 1), U a (p, 4) = u a (p, 2). (A3) 
u and v are given as follows with Pauli matrix a and two dimensional basis vector £(r), 



v a (p,r) 
u a (p,r) 



\X'(p)\-B(p) 



\ 2|A'(P)| \w^M\^r) 



a 



1, -io 

1 



\ 2|A'(p)| 
, ^(2) = 



(A4) 

(A5) 

(A6) 
(A7) 
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U a (p, s) satisfies the orthogonal and completeness condition 

"ZUl{p,s)U a (p,t) = 6 Stt , (A8) 

a 

J2ut(p,s)U fj (p,s) = 5 a ,p. (A9) 



APPENDIX B: ANOMALOUS QUARK MASS WITH COMPLETELY 
LOCALIZED AND PLANE- WAVE EIGENVECTORS 

We consider a system with jVj completely localized eigenvectors and plane- wave functions 
with N p degrees of freedom for momentum p. 

We expand the two Green functions in m 5g with the complete set in eq.( [4.35| ). 

X(t) = J2 (J 5q (x)P(y)) 



Q5 

n 10 



EE 7 



x,a,a b,f3 



2 cosh ^a 5 H 



J) J 



K 1 )' 



2 cosh ^a 5 H 



«5 



x ( J 



( Ni 

EE E m 

x,a,a b,/3 \n,m=l 



n 



2 cosh 4^a 5 # 



+ E E ( J b> s > c ) (p> s , c 

P,»,ck,t,d \ 



2 cosh ^a 5 i/ * 8 
(m|J) 

1 j 



J 



x ( J 



2 cosh ^a 5 # 



2 cosh ^a 5 H Ns 



k,t,d)(k,t,d\I) 



where / = (x, a, a), J = (y, (3, b) and we use a relation (/|n) (A;, t, d\I) = 0. 
Contribution from the localized eigenvector is given as follows 



Ut) = -^EEE (m Hi) 



x,a,a b,P n,m=l 
1 



x n 



2 cosh ^a 5 H 



J) J 



D ~nI 



2 cosh ^-a§H 



rn 



Ni 



-^0 E 5 *,^ 



n=l 



v E( n d nI y,f3,b) {y,P,b 



(Bl) 



n 



where we use a relation 



(I\n) = (x,a,a\n} = 5 a>an 5 a ,a n 5 x 
(I\n) (m\I) = S mn 8 aan 5 aan 5 XXn 



(B2) 

(B3) 
(B4) 
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By inserting the identity operator in (x, a, a) space 

1 = \x, a, a) ( 



(B5) 



we have 



a 5 



Ni 



a 



10 



E s t,x° 



2cosh^a 5 A(A0 



x E E E ( n k>7,c; < :. -.c 

b,(3 z,7,ciu,5,d 



y,P,b)(y,P,b 



w,5,d) (w, 8, d\n) 



a 5 ^ 
-Zva E 



1 



a 



n=l 



2cosh^a 5 A(A{J 



n 1 N < 

'"" E 



1 



a w AN C N X ^ ^coshf a 5 A(A0 
x E E(^ a 'H^1 y '^' & ) 

2 



1 



= -Y(t) 



E 



a 5 47V C 7V X ^ I 2 cosh f a 5 A(A; t ) i ' 



(B6) 



where we assume that the center of the localized solution (x n ,a n ,a n ) is distributed uni- 
formly in the (x, a, a) space and we can take average over (x, a, a) when enough number of 
configurations are summed in the simulation. 

For the contribution from the continuous eigenvalues we take t — > oo limit to proceed, 



a x,a,a b,/3 P,s k,s> ly P 



x(p,s,a 



2 cosh ^a 5 H 



--DZI 



y,P,b)(y,P,b 



K) 1 



2 cosh ^a 5 H 



k, s', a 



a 5 i J_ \ - 

~^ n ^^2coshfa 5 A(p) 
xJ2(p,s,a\DZ 1 5 \y,P,b)(y,P,b (d^) 1 



p,s,a), 



where we use following relations 

(x, a, a\p, s, b) 



-U a (p,s)e^8 atb , 



? 6 

X 



.(p-k)x _ I r _ 



(B7) 

(B8) 
(B9) 
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lim e lt{p °- k0) = <Lo fc0 , 



(BIO) 
(Bll) 



is a volume of 5 in xyz space. By inserting the complete set ( |4.35|) of eigenstates we 

have 



*e(*)| 



t— >oo 



a 10 xyz N pp %\ 2 cosh fa 5 A(p) 



E 



X E( E \^ S ' a 



\n,m=l 

+ E (p 



s, a 



( n \y,f3,b) (y,(3 } b\m) (m 
h,s',c) (k } s',c\y,{3 7 b) 



p, s, a 



X E {y^M^s\d) (q,s",d 

g,s",d 



p, s, a 



Q5 , J_ 



E 



1 



2 cosh -^az,\(p) 



x 



E E <Wn \P> S > 
\ a,s n =l 



n) (n 



K 1 )' 



p, s, a 



05 / J_ \ - 

■^ n T„Z N 2^ 



C 2 + E 2 
1 



/ » 7fi C M (p) + ^(p) \ 



s,s' 
2 



P p 



2cosh^a 5 A(p) 



F(p), 



where F(p) is defined as 



F(p) = F l (p) + F c (p), 



Ni 



F 1(P) = E <W» E ^2 \ P ' S ' ° 



n=l 

4iV, 



n) (n 



p, s, a 



N p \C 2 (p) + E 2 (p) J ' 



(B12) 

(B13) 
(B14) 

(B15) 



We can show that Fi(p) vanishes as follows. We extend the orthogonality assumption of the 
eigenvectors (n\p,s,a) = to the case with 75; (n\j5\p,s,a) = 0. This is plausible since 
the eigenfunction ip n (x) = (x\n) is localized in space-time and the overlap is suppressed 
even if we multiply 75. Then by using the explicit form ( |3.5|) , the matrix element of the 
truncated overlap Dirac operator D^ 5 turns out to be block diagonal, in which only the 
matrix elements (n \D N5 \ m) and (p, s, a \D Ns \ k, t, b) are non-zero and the off-diagonal part 
{p, s, a \Dn 5 \ n) becomes zero. This block diagonal property is kept even if we take inversion 
and the off-diagonal part becomes zero 



p, s, a 



DnI 



?? 



0. 



(B16) 
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It is noted that in the free theory the relation D^Dw = DwDw is satisfied and 



n 



as 



D N * mj = is shown exactly for non-degenerate \X' n \ ^ \X' m \. 

The pion propagator in the denominator of m 5g is expanded similarly with eigenstates 



Y(t) = £ (P(x)P(y)) 



,10 



EE( J d- n i 



x,a,a b,/3 



J) J 



|EE(E (I\n)(n\D- N l 

x,a,a b,/3 \n,m=l 



P,s,ck,s',d 



(D^)t m) <m|/> 



(£>^)t fc, s ',d)(fc, s / ,d|/> 



(B17) 



A contribution from localized eigenvector is given by 

of * 
a 1 



Y iit) = - E ^ E (« ^ V> ft 6) ft b TOj)t 
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= — fo E 5 M« E E H z > 7, c) (2, 7, c AvJ y, ft & 
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w,S,d 



a 2 Nl 

-^E^E(^ 
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4N C N X 



y,/3,b)(y,/3,b (D^Y 
E ^(x,t,a,a\D]£\y,P,b)(y,P,b\(Djfiy x,t,a,a ; 

c,a,a (3,b 

Y(t). (B18) 



io ^ AN C N X ^ jri 



Taking t — > oo limit a contribution from the continuous sector becomes 
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+ E \Pi s i a d n 5 k,s',c) (k,s',c\y,/3,b) 



k,s'c 



x £ (y,(3,b\q, S ",d)(q, S ",d (D^ 



p,s,a 



(B19) 



As a consequence we have a relation 



Y{t -> oo) 



a| 1 



4AU\L - iV, a 10 ^iV, 



(B20) 



p p 



Now we need to know behavior of F(p) as a function of p to further proceed. We can 
directly see that the function F(p) = F c (p) is almost constant for all the range of p as in 
Fig. 0. We assume that F(p) is independent of p. 

By adopting this assumption the anomalous quark mass becomes 



m 5q 



lim^ 00 F(t) 

1 1 E 



a 5 4:N C N X „ I 2 cosh 4fa 5 A(A^ 



° 5 ' J_ \ "* 
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E 



a 5 4iV,JV x V \ 2 cosh f a 5 A(A; 

1 



AN C N X - Ni 1 



E 



AN C N X a 5 E fc V I 2 cosh f a 5 A(p) 



F(p) 
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E 



a 5 4N c N x x „ 
1 1 



2cosh^a 5 A(Ag 



a 5 ±N C N X 



y 

^ V 2cosh^a 5 A(A'„ 



+ Ep« 



2cosh^a 5 A(jo) 
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(B21) 



where p n is a number of degeneracy in A n , 

Pn = Yl * 



A n ,A(p)' 



(B22) 
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APPENDIX C: ANOMALOUS QUARK MASS WITH PARTIALLY LOCALIZED 

AND PLANE- WAVE EIGENVECTORS 



We consider a system with Ni subspaces S n (n = 1, • • • , iVj), which is spanned by partially 
localized eigenvectors ip$. We assume that each subspace S n has no overlap and its volume 
is h n . Remaining subspace S of the system is spanned by plane- wave eigenvectors as in the 
previous subsection. If we set the size of S to be 4N c N p the total degrees of freedom of the 
system becomes 



AN C N X = J2 h n + 4N C N P . 



n=l 



The partially localized eigenvector is given by 

H'\n,i) = \W\n,i), 
(I\n,i) = 4=V£°(J), 



'hr, 



where the function ipw(I) is non-zero within S n . We have a relation that 

#(i)tf(/) = 
for n ^ m for any The orthogonality is given by 



(CI) 



(C2) 
(C3) 



(C4) 



(C5) 



The plane- wave eigenvector is the same as in the previous section ( 4.25| ) with a constraint 



<n,i|/> (l\p,s,b) = 0. 
The completeness condition of this system is given by 



n=l i=l 



p a 



with 



Ni hn I 

h, 



n=l i=l 

We expand the two Green functions in m^ q with the complete set in the above. 

Mt) = H (Mx)P{y)) 



4EE ' 
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2 cosh ^a 5 # 



j j 
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2 cosh ^a 5 H 
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+ E E ( x ' a > °b> s > c ) (p> s > c 
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2 cosh ^a 5 i/ ' ' 



J 



x J 



2 cosh ^a 5 # 
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The contribution X c from the continuous eigenvalues is the same as in the previous subsection 
except for the definition of F(p) 



F(p) = F l (p) + F c (p), 

Ni h n i / , 

Flip) = E E E(p.MNh.i>r#Wf ) W(n.J (^;) 
47V / 4 \ 



p,s,a 



N p \C 2 (p) + E 2 (p) J ' 
Here we use the block diagonal condition 



p,s,a 



DnI 



n,i) = 



(CIO) 
(Cll) 

(C12) 
(C13) 



by assuming the 75 orthogonality (n, i \ j 5 \ p, s, a) = 0. X c is written in terms of the function 
F{p) = F c {p). 

Contribution from the localized mode becomes 

^)-4EEEr/: j » 

u b,/3 n=l ij=l lbn 



X (n,t 



2cosh^a 5 # ^ 



J J 



2 cosh ^a 5 # 



n Ni h n 1 

u n=l jj=l n n 



1 



2cosh^a 5 A(A^ ) ; 



2 cosh f a 5 A(An 0) ; 



x E E E ( n ' 7, c) 7, c D^J 

„ Ni h n -t 

a.5 f-^ 1 



J) ( J 



u n=l jj=l 
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2cosh^a 5 A(Al W ; 



2 cosh f a 5 A(An 0) ; 



x E E E r-^ 0t (*> 7, c) (z, 7, c |l>^ 



6,/3 2,7^111,(5,(1 ^™ 



J J 



K) f w,6,d)1,S\w,6,d), 

(C14) 

where we assumed orthogonality for m and n, is defined as 
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■,(j)t 



(C15) 



x a, a 



fn^{t) is non-zero when t e «S n . Here we have a comment. It is plausible to assume that the 
eigenfunction ip^ tends to be plane but rapidly oscillating inside S n for large i. By using 
this fact we can see that the above is suppressed for different i, j even if the summation 
over t is not taken. 

We need more information on the propagator multiplied with eigenvectors 



1 



n I K 



J) J 



(d^\k)^\K) 



(C16) 



to further discuss the detailed property of X[. We start by adopting an assumption that the 
propagator 



J) J 



K 



(C17) 



is a slowly varying function of I and K. Then we investigate the behavior of (|U16|) for two 
typical forms of eigenf unctions. The eigenfunctions are classified into two types, according 
to the behavior of the following integral of a single function 



c = £V#(/)- 



(C18) 



(i) Eigenfunction ipffi with C ^ 0, which may be typical for the lowest mode i = 1. It 
seems to be plausible |T9|J3T| to approximate this lowest mode with the exponential form. 
For simplicity we consider a one- dimensional case 



if>(x) = — e -\ x - X0 U 5 , 



(C19) 



where 5 is a width and xo is a center of the eigenfunction. Since the integral / dxip{x) gives 
non-zero value, a multiplication with a smooth function f(x,y) produces a enhance factor 
45, which is proportional to a width of the eigenfunction: 



1 



dxdye- lx - Xol/5 f(x, y)e~ ly - Xol/5 



1 r°° 

- / dzdwe~ lzl/6 f(z + x ,w + 

J-oo 



x )e 



— | ui I / 8 



-oo 

6 



dzdwe-^ 5 f(x ,x )e-^ s + O(df(x ,x )) = (45) f(x ,x ). 



(C20) 



The extension to four-dimensions is given by 



^\x) 



J2i \ x i~( x n)i\/ S 



(C21) 



and (|C16|) becomes 
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ll"n j k 
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where I n is a peak of the exponential in Vri (-0 an< ^ ^« ^ s gi yen by fo n = <5 4 - Here we use a 
property that the eigenfunction ip^ with higher i > 1 tends to be oscillating and a single 
summation Y<i (-0 is suppressed. Together with the suppression of /* J for i ^ j we have 
a factor 5^ in the above. 

The term X\ of the numerator for the contribution from the lowest eigenvector i = 1 
becomes 



n=l i=l 



2cosh^a 5 A(A„ (i 
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a 5 4Ar c iy x ^ \ 2 cosh ^a 5 A(A 



(C23) 



Here we use a relation for i = j 



7 " ^tU n ? 



(C24) 



where /i n) i is a width of t in «S n . 

(ii) Plane but rapidly oscillating eigenvector with C ~ 0, which is typical for the 
excited mode with i > 1. In this case as was mentioned in the above a single summation 
X)/ (I) is suppressed and ( |C16| ) has non-zero value only when I = K, 



f ££#00(1 
1 £# f (/)(/ 



/?., 



IW<P(I) 
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where /, J unsummed. Xi then becomes 
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In this case m 5q is given by the same formula as in the previous section with the completely 
localized eigenvector. In addition we assume that the spin and color eigenfunction is always 
rapidly oscillating and set their enhance factor to be unity as is discussed in the above. 
The pion propagator in the denominator is expanded as 

Y(t) = £ (P(x)P(y)) 



„2 / N[ h n -l 

4EE E E^ l) W^(/)(v 

u x,a,a b,/3 \n,m=l i,j=l n n 



DnI 
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>,d 1S P 
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By inserting a complete set 1 = 1-0 (I\ the contribution from the localized eigenvectors is 
rewritten in term of Y(t). 
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n 2 _ 111 hn 1 / 
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(C28) 



where we take average for I e U n iS n . 

The contribution from the plane-wave modes is the same as that of the completely 



localized case by using the off diagonal property (p, s, a 
relation to write down Y(t) in terms of F(p) only, 



D N 5 



n, i 



0. Now we have a 



Y(t -> oo) = 



AN r N x 



a\ , 4iV c 



47V c iV x - hty a 10 *»* iV p ^ 



(C29) 



where /iA^ = Y^^Li h n , and the anomalous quark mass m 5q can be written as 
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lim f ^ 00 X(t) 



Ni ( 1 \ 2 N t h n 

... r„„J + ££ 



" ^2 cosh f 05 A(A; (1 M „t- 1 Kl2coshfa 6 A(A: (,) ) 



2. 



+ c ^2 Pn (^2cosh^a 5 A n J j ^ ^ 

In the actual case, the weight factor for the contribution from partially localized eigen- 
states lies between 1 and A A h n depending on the index i. Therefore we adopt h$ such that 
1 < h$ < A 4 h n as the weight factor in the text. 
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FIGURES 
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FIG. 1. The expected phase structure of lattice QCD with the Wilson fermion in the M-g 2 

plane. In the region A, (ijji'y^T^tp) = (the parity- flavor symmetric phase), while (vpi'jsT 5 ^} ^ 

(the parity- flavor broken phase) in the region B. 
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FIG. 2. F c (p) as a function of A(p) at M = 0.8(left) and M = 1.4(right). Both F c and A 
are evaluated on 16 3 x 32 momentum lattice, where the periodic boundary condition is adopted 
for spatial direction and anti-periodic boundary condition is set for temporal direction in order to 
prevent a singularity at the origin. 
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FIG. 3. Three types of p(X) adopted in calculation of m^ q . The left one corresponds to type 
(1) with 5 < 0, the center is type (2) with 5 = and the right is type (3) with 5 > 0. 
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FIG. 4. m,5q as function of N$ evaluated with three types of p(X). The left one corresponds to 
type (1) with 5 = —0.5, the center is type (2) with 5 = and the right is type (3) with 5 = 0.5. 
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